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Abstract. This is the fourth paper of our study of the shape parameter c contained in the fa- 
mous muhiquadrics (— l)^'^! (c-^ + /3 > 0, and the inverse muhiquadrics (c^ + (3 < 0. 
The theoretical ground is the same as that of [TU]. However we extend the space of interpolated 
functions to a more general one. This leads to a totally different set of criteria of choosing c. 

keywords: radial basis function, multiquadric, shape parameter, interpolation. 

1 Introduction 

Again, we are going to adopt the radial function 

h{x) := r(-^)(c2 + \x\'')^, 13 € i?\27V>o, O (1) 

, where |a:;| is the Euclidean norm oix in i?", T is the classical gamma function, and c, /3 are constants. 
This definition looks more complicated than the ones mentioned in the abstract. However it will 
simplify the Fourier transform of h and our analysis of some useful results. 

In order to make this paper more readable, we review some basic ingredients mentioned in the 
previous papers, at the cost of wasting a few pages. 

For any interpolated function /, our interpolating function will be of the form 

N 

s(x) -.^^^a^x - Xi) + p{x) (2) 
1=1 

where p{x) G Pm-i, the space of polynomials of degree less than or equal to m — 1 in i?", A — 
{xi, ■ ■ ■ ,xn} is the set of centers(interpolation points). For m = 0, F„i-i '■= {0}. We require that 
s(-) interpolate /(•) at data points (xi, /(xi)), • • • , (xjv, /(a^iv))- This results in a linear system of 
the form 

N Q 
^C^h{xJ - X^) + ^b^p^{xJ) = f{Xj) , j = 1, • • • , A 



(3) 
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N 

^CiPj(x^) ^0 ,j = l,...,Q 

to be solved, where {pi, ■ ■ • ,pq} is a basis of Pm-i- 

This hnear system is solvable because h{x) is conditionally positive definite (c. p. d.) of order 
m = max{[^],0} where \^~\} denotes the smallest integer greater than or equal to ^. 

Besides the linear system, another important object is the function space. Each function of the 
form (1) induces a function space called native space denoted by Ch,m{R"), abbreviated as Ch,m, 
where m denotes its order of conditional positive definiteness. For each member / of Ch,m there is 
a seminorm called the /i-norm of /. The definition and characterization of the native space 

can be found in [4] , [5] , [7] , [11] , [12] and [14] . In this paper all interpolated functions belong to the 
native space. 

Although our interpolated functions are defined in the entire i?", interpolation will occur in a 
simplex. The definition of simplex can be found in |3]. A 1-simplex is a line segment, a 2-simplex is 
a triangle, and a 3-simplex is a tetrahedron with four vertices. 

Let T„ be an n-simplex in i?" and v^, 1 < i < n + 1 be its vertices. Then any point x G T„ can 
be written as convex combination of the vertices: 

n+l n+1 



= ^ CiVi, ^ Cj = 1, Ci> 0. 



i=l i=l 

The numbers ci, • • • ,c„4.i are called the barycentric coordinates of x. For any n-simplex T„, the 
evenly spaced points of degree / are those points whose barycentric coordinates are of the form 

(-^, • • • , ~y~")i nonnegative integers with ki = I. 

i—l 

It's easily seen that the number of evenly spaced points of degree / in T„ is exactly 

N = dimPP = f " + ^ 
' \^ n 

where P" denotes the space of polynomials of degree not exceeding I in n variables. Moreover, such 
points form a determining set for P", as is shown in [5]. 

In this paper the evaluation argument x will be a point in an n-simplex, and the set X of centers 
will be the evenly spaced points in that n-simplex. 

2 Fundamental Theory 

Before introducing the main theorem, we need to define two constants. 
Definition 2.1 Let n and 13 he as in (1). The numbers p and Ag are defined as follows, 
(a) Suppose /3 < rt — 3. Let s = [ "~^~'^ ] . Then 

(^) zf(3<0, p=^^andAo = ^2±M^±^; 
(U) z//3>0, p=l + ^ and Ao = (2»^+^+^)(2;»+j+^)-(2»»+3) 
where m = [^] . 
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(b) Suppose 71 — 3</3<n — 1. Then p = I and Aq = 1. 



(c) Suppose /3 > n — 1. Let s = — [ " ^ ^ ] . T/ien 

p = 1 and Aq — r where m = [— 1 . 

^ " (2m + 2)(2TO+l)---(2?7i-s + 3) '2' 

The following theorem is the cornerstone of our theory. We cite it directly from |5] with a slight 
modification to make it easier to understand. 

Theorem 2.2 Let h be as in (1). For any positive number bo, let C ~ max | gf^, 8p^ and Sq = 
For any n-simplex Q of diameter r satisfying < r < (note that ^ < bo), if f & C/i.„j, 

|/(.t) - s{x)\ < 2'^7T'^^TI^,c^-'^V3CVS{y)^\f\\h (4) 

holds for all X ^ Q and Q < 5 < 5q, where s{x) is defined as in (2) with xi, • • • , xjsi the evenly spaced 
points of degree I in Q satisfying < I < . The constant an denotes the volume of the unit 
ball in R", and < A' < 1 is given by 

which only in some cases mildly depends on the dimension n. 

Remark: (a)Note that the right-hand side of (4) approaches zero as (5 — )• O"*". This is the key to 
understanding Theorem2.2. The number 5 is in spirit equivalent to the well-known fill-distance. 
Although the centers Xi, - ■ ■ ,xn are not purely scattered, the shape of the simplex is controlled by 
us. Hence the distribution of the centers is practically quite flexible. (b)In (4) the shape parameter 
c plays a crucial role and greatly influences the error bound. This provides us with a theoretical 
ground of choosing the optimal c. However we need further work before presenting useful criteria. 
In this paper all interpolated functions belong to a kind of space deflned as follows. 

Definition 2.3 For any positive number a, 



:= l^f e L\Rn : J \f{0\'e^ d^ < 



where f denotes the Fourier transform of f . For each f G E^, its norm is 

1/2 



The following lemma is cited from [9]. 

Lemma 2.4 Let h be as in (1). For any u > 0, i/ /3 < 0, |n + /3| > 1 and n + P + 1 > {), then 
Ea ^ Ch,m{R^) and for any f G E^;, the seminorm \\f\\h of f satisfies 



whe 



1/2 



ccr + y/c'^a^ + 4cr(n + 13 + 1) 
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Corollary 2.5 Under the conditions of Theorem2.2, if f £ E^, /? < 0, |n+/3| > 1 andn + fi + l > 0, 
(4-) can be transformed into 

\f{x) - s{x)\ <2~'^-^7:-i--iy^y^oV3Cc^^^^ i^ia^e'^'-^^^^^ (5) 
where 

_ ca + ^cV^ + 4a(n + 13 + 1) 



Proof. This is an immediate result of Theorem2. 2 and Lenima2. 4. jj 

Note that Corollary2.5 covers the very useful case /3 = — 1, n> 2. However the case /3 = — 1, n = 1 
is excluded. For this case we need a different approach. 

Lemma 2.6 Let a > 0, /3 = —1 and n = 1. For any f e Ea-, 

1/2 



where M{c) := e^"^ if c < andM{c) := g(££+v^i!+ii) if c > t/;/iere g(^) := ^/^e= 



p^S 

/3o" 

Proof. This is just Theorem2.5 of [S]. jJ 

Corollary 2.7 Let cr > 0, /3 = — 1 and n= 1. Under the conditions of Theorem2.2, if f £ Ea-, (4) 
can he transformed into 

1/2 



\f{x) ~ s{x)\ < ^^^VsCc^-' I — + 2V3M{c) [ V^(A')*|1/|| 



(6) 



where M(c) is defined as in Lemma2.6. 

Proof. This is an immediate result of Theorem2.2 and Lcmma2.6. 

Now we have dealt with the most useful cases for (3 < 0. The next step is to treat /3 > 0. 
Lemma 2.8 Let a > 0, (3 > and n> 1. For any f £ Ea, 



1/2 



u<rfoc-^ r- ' ^'^^ I \\f\\E, 



where ^* = '^"'+V'^ °" +4g-(i+ff+") ^ constant depending on n, (3 only. 

Proof This is just Theorem2.8 of [3]. 

Corollary 2.9 Lei ct > 0, (3 > and n> I. If f £ Ea, (4) can be transformed into 



\f{x) - s{x)\ < 2^7T'^y/E^^oV3Cdoc'^^^ I , ' , |> V^(A')*||/b. (7) 



1/2 



e 



where do, ^* are as in Lemma2.8. 

Proof. This is an immediate result of Theorem2.2 and Lemma2.8. 
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3 Criteria of Choosing c 

Note that in (5), (6) and (7), there is a main function of c. As in [9 , let's call this function the 
MN function, denoted by MN{c), and its graph the MN curve. The optimal choice of c is then the 
number minimizing MN{c). However, unlike (H], the range of c is the entire interval (0, oo), rather 
than a proper subset of (0, oo). 
We now begin our criteria. 



Casel. 



P <0, |n + ;3| > 1 and n + /3 + 1 > Let f £ Ea and h be as in (1). Under the con- 
ditions of Theorem2.2, for any fixed 5 satisfying < (5 < Jq, the optimal value of c in (0, oo) is the 
number minimizing 

1/2 



ff-„ + l-4i I n + g + 1 ({*) 



MN{c) := c J <^ iC) 



e 



where 



+ ^c2cr2 + 4cr(„ + ^ + 1) 



Reason: This is a direct consequence of (5). j) 

Remark: (a)It's easily seen that MN{c) — oo as c — >■ oo. Also, if n + /3 + 1 > 0, MN{c) — )■ oo 
as c — > 0+. (b)Casel covers the frequently seen case /3 = —1, n > 2. (c)The number c minimizing 
MN{c) can be easily found by Mathematica or Matlab. 

Numerical Results: 



Graph of the MN Curve with 5=0.01 and 1=5 

MN(c) 
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Graph of the MN Curve with 5=0.008 and 1=6 

MN(c) 




Figure 2: Here n = 2, /3 = —1, a ~ 1 and = 1. 
Graph of the MN Curve with 5=0.006 and 1=7 

MN(c) 



0.0015 - 



0.0010 - 




Figure 3: Here n = 2, (3 — —1, cr = 1 and bo = 1. 
Graph of the MN Curve with 5=0.004 and 1= 1 1 

MN(r) 




Figure 4: Here n = 2, l3 — —l,a = I and bo ~ 1. 
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Graph of the MN Curve with 5=0.002 and 1=21 

MN(c) 
1.x 10-" - 



8.x 10-"* 



6.X10- 



4. X 10-"* - 




Figure 5: Here n = 2,(3 — —1, a — I and bo — 1. 



Case2. /3 ~ —1 and n — 1 Let f £ and h be as in (1). Under the conditions of Theorem2.2, 



for any fixed 5 satisfying < 5 < (5o, the optimal value of c in (0, oo) is the number minimizing 

1/2 



MNic) := c^-' I A -I- 2V3M(c)l 
in2 J 



where 



M{c) :-- 



1 — i 



if < c< 



„i- ccr+Vc^g-^+4g- ^ -r 2 ^ 



, g being defined by g{£) := ^/c^e'^^~~ . 

Reason: This is a direct result of (6). jj 

Remark: Note that MN{c) — > oo both as c — ?• oo and c — t- 0+. Now let's see some numerical 
examples. 

Graph of the MN Curve with (5=0.01 and 1=5 

MN(c) 




Figure 6: Here n = l,/3 = — l,cr=l and 6o = 1- 
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Graph of the MN Curve with (5=0.008 and 1=6 

MN(c) 




4 6 8 10 



Figure 7: Here n = 1, l3 = —l,a — I and bo = 1. 
Graph of the MN Curve with 5=0.006 and 1=7 

MN(c) 




4 6 8 10 



Figure 8: Here n = l,/3 = — 1,<t=1 and bo — 1. 
Graph of the MN Curve with 5=0.004 and 1= 1 1 

MN(f) 




Figure 9: Here n = 1, (3 — —l,a ~ 1 and bo = 1- 



8 



Graph of the MN Curve with 5=0.002 and 1=21 

MN(c) 



3.x 10-" 
2.5x10-" 

2.x 10-" 
1.5x10-" 

1.x 10-" 
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12 



Figure 10: Here n — l,/3~ — l,cr = 1 and bo = 1. 



Cases. 



(3 > and n > 1 Let f £ and ft, be as in (1). Under the conditions of Theorem2.2, 



for any fixed 5 satisfying < (5 < (5o, the optimal value of c in (0, oo) is the number minimizing 



MN{c) 



i+g-,i-4i I (4 ) 2 I 



e " 



, where 



CO- + y^c2cr2 + 4cr(l + /J + 



Reason: This follows from (7) 
Remark: By observing that 



1 r 



16 



2cV + 2cVc2cr2 + 4cr(n + /? + 1) - (4?i + /3 + 1) 



, we can easily obtain useful results as follows. (a)If 1 + /3 — n — 4Z>0, limc_i.o+ MN{c) = 0. (b)If 
l + P-n-4:l<0, limc_>o+ ^(c) = oo. (c)If l + /3-n-4Z = 0, hmc^o+ M N{c) is a finite positive 
number. (d)limc_>oo MN{c) = oo. 

Numerical Results: For simplicity, we offer results for n = 1 only. In fact for n > 1 similar 
results can be presented without slight difficulty. 
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Graph of the MN Curve with (5=0.01 and 1=5 

MN(c) 
2.0 - 



1.5 - 
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Figure 11: Here n=l,/3=l,cr = l and bo — 1. 
Graph of the MN Curve with (5=0.008 and 1=6 

MN(c) 
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Figure 12: Here n = 1, f3 = l,a = 1 and 6o = 1- 
Graph of the MN Curve with 5=0.006 and 1=7 

MN(c) 



0.08 



O.P6 - 



0.04 - 



0.02\ 
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Figure 13: Here n = 1,^= 1, cr = 1 and bo ~ 1. 



Graph of the MN Curve with 5=0.004 and 1= 1 1 

MN(c) 




Figure 14: Here n — 1, l3 = l,a = I and bg = 1. 
Graph of the MN Curve with (5=0.002 and 1=21 

MN(c) 




Figure 15: Here n = l,/3=l,(T = l and bo ~ 1. 
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